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Abstract 

We derive a universal function for the kinetics of complex systems characterized by 
stretched exponential and/or power-law behaviors. This kinetic function unifies and 
generalizes previous theoretical attempts to describe what has been called "fractal 
kinetic" . 

The concentration evolutionary equation is formally similar to the relaxation func- 
tion obtained in the stochastic theory of relaxation, with two exponents a and n. 
The first one is due to memory effects and short-range correlations and the second 
one finds its origin in the long-range correlations and geometrical frustrations which 
give rise to ageing behavior. These effects can be formally handled by introducing 
adequate probability distributions for the rate coefficient. We show that the distri- 
bution of rate coefficients is the consequence of local variations of the free energy 
(energy landscape) appearing in the exponent of the Arrhenius formula. 

The fractal (n, a) kinetic is the applied to a few problems of fundamental and 
practical importance in particular the sorption of dissolved contaminants in liquid 
phase. Contrary to the usual practice in that field, we found that the exponent a, 
which is implicitly equal to 1 in the traditional analysis of kinetic data in terms of 
first or -second order reactions, is a relevant and useful parameter to characterize the 
kinetics of complex systems. It is formally related to the system energy landscape 
which depends on physical, chemical and biological internal and external factors. 

We discuss briefly the relation of the (n, a) kinetic formalism with the Tsallis 
theory of nonextensive systems. 
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Introduction 

The physical origin of anomalous kinetics in complex systems like glass, poly- 
mers, solutions, proteins, enzymes... has received much attention in recent 
years ([1][2][3][4][5]). 

This is due to the observation that in many instances, the kinetics cannot 
be accounted for without introducing a time-dependent reaction rate coeffi- 
cient to describe properly the experimental data. As in many other related 
problems an accurate relation between microscopic properties and the global 
macroscopic observables is lacking due to many-body interactions, the in- 
evitable homogenization and coarse-graining resulting from complexity and 
experimental techniques. What has been done in practice to deal with this 
situation is to generalize formula used in simple reaction kinetics (first and 
second order reaction) and introduce one or several supplementary empirical 
parameters to fit experimental data. The aim of these works is to establish 
correlations between macroscopic observables and external parameters (tem- 
perature, concentration, pH, ...). Another method widely used is to consider a 
few reaction steps. The total rate equation is written as the sum of elementary 
rate equations (first or second order) . Their respective weight and the rate co- 
efficients are fitted to the data. This procedure can lead to ambiguous results 
and conclusions. In particular it does not reproduce power-law behaviors often 
encountered in these systems. In this spirit, Frauenfelder and collaborators [3] 
have used intuitive arguments to fit experimental biomolecular reaction data 
in protein materials with empirical stretched exponential or power-law func- 
tions. They trace the 'anomalous" kinetics to the distribution of the reaction 
rate in the Arrehnius formula. In 1988, to account for power law behavior, 
Kopclman [1] proposed a phcnomenological fractal like kinetics to account for 
reactions in materials prepared as fractals. This lead more recently Savageau 
[2] to introduce a model where instead of introducing a time-dependence to the 
rate coefficient, the reactant concentrations are raised to non-integer powers. 

More recently, Weron [4] and collaborators using results from their stochas- 
tic theory of relaxation in dipolar systems [6], introduced a time dependent 
power-law reaction rate coefficient to generalize first and second order kinetic 
equations in order to apply them to biomolecular reactions. Simultaneously, 
Mendes and collaborators [7], using results of the Tsallis nonextensive entropy 
theory [8] to solve non-linear differential equations, introduced the concept of 
a n-order kinetic equation whose solution has a formal expression similar to 
the TsaUis generafized Pareto distribution. 

The purpose of the present paper is to use some results of two recent papers on 
non-Debyc relaxation [9] [10] to incorporate the ideas developed in the previous 
quoted works [1][3][4][7], in one unified formalism in order to the set the basis 
of a general theory of reaction kinetics in complex systems. 
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The challenge which has to be faced, in this important problem, is to give a 
physical or statistical meaning to the empirical parameters and, when this is 
possible, to relate the "anomalous" kinetics to universality, i.e. universal scal- 
ing functions independent on the microscopic or mesoscopic detailed properties 
of the system. 



1 The (n, a) kinetic equation 



The most general chemical kinetic equation for one given species (^4) in a 
complex system composed of A,B,C,.. reacting atoms and molecules can be 
written formally as 

^ = Kcldcl... (1) 

where K is the rate coefficient and a, /3, 7. ..refer to the concentrations of chem- 
ical species A, B, C, ...present in the reaction and the sum a + f3 + 'j is the 
overall order of the reaction. In some cases, the concentration CbCc--- can be 
considered as constant, thus the above equation reduces (for instance for the 
reactant A) to the form 

^'"^^^ - K^c^itr (2) 



dt 

In that way, the parameter n becomes the overall order of the reaction. The 
solution of this differential equation for one of the reactants is given by [7] 

Cn{t) = c„(0)[(l + (n - l)c„(0)"-^X„t]-^ (3) 

which has the form of a generalized Pareto function, solution of the Tsallis 
entropy maximization [8] and has an asymptotic power law behavior c„(t) oc 

If we use the deformed n-exponential and n-logarithm introduced by Tsallis 
and collaborators in the context of nonextensive systems [8]: 

exp„(x) — {1 — {n — if 1 — {n — l)x > 0, otherwise 



and ln„(x) = — (4) 

1 — n 



with 

exp„( hin{x)) =lnn( exp„(a;)) = 1 
we can write the solution (3) in a more compact form: 



Cn{t) = c„(0) exp„(-t/rn) (5) 
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with a characteristic time: 

Tn = {c^m'-^K^r^ (6) 

This will allow the definition of an effective time-dependent rate coefficient: 
^^^'^ = c(0)l(expJ-(Vr„))" (7) 



dt T,_ 

= c(0)-(l + (n-l)(t/r„))-^-^ 

= lCn{t)Cn{t) 

with 

/C„(t) = -(l + (n-l)(VT„))-^ (8) 
For t «Tn, one has a slowing down of the effective rate : 

/C„(i) = -(l-(n-l)(t/r„)) t«T„ (9) 

and for n 7^ 1, t >> t„ 

/C„(t) (X (l/(n - l))r^ (10) 

This behavior is a manifestation of what has been call aging [11] which appears 
as soon as n 7^ 1. For n — 1, one recovers the exponential behavior with 
lC{t) = l/r. 

These results do not exhibit the t <t power law time dependence of the reac- 
tion rate which describes adequately the experimental data of many complex 
systems [4] [12] [13] [14] [15]. This behavior can appear quite naturally if we in- 
troduce in (6), instead of the n-exponential, a n— Weibull function introduced 
by Mendes [16] and used also in the theory of relaxation [9] [10] 

Cn,a{t) = C(0) expJ-(t/T„,„)") = c(0)[(l + (n - 1) (VTn,a)"]~^ (H) 

with a characteristic time: 

Tn,a = [Xn,aC(0)'^-^]-^/« (12) 

The effective time-dependent rate coefficient Kn,aii) now reads 

^n,a{t) = + - m/rn,ar)-' (13) 

Equation (11) is solution of a fractional differential equation : 
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by introducing a fractional time index a and a non-integer reaction order 
n. Fractional derivation and fractal time concepts have been introduced in 
physics (diffusion in disordered and porous media, random walks •••[17]) in 
the theory of dielectric response [20] [18] [10] and in economy [19] • 

The effective rate coefficient }Cn,a{t) has the two asymptotic behaviors 

for t^O }Cn,a{t) « ^""^ (15) 

for i — > OO K.n,a{t) oc t"^ 

For t ^ 0, we get the same power-law variation of the rate coefficient as in the 
work of Weron et al. [4] as well as in the fractal phenomenological description 
of non-homogeneous reaction dynamics called fractal-hke kinetics [1], if we 
identify the Kopelman fractal parameter h < 1 with 1 — a. As noted in [5], 
the concept of effective time-dependent rate constant breaks down for t — > 0, 
a < 1, since in that limit /C„,q,(^) diverges. The general solution (11) of the 
fractional differential equation (14) does not suffer from such difficulty and is 
well defined in the positive time domain. In any case, as for real geometric 
fractals, for physical reasons, there is in each case a natural small time cut-off. 

The two asymptotic behaviors of the concentration evolutionary law equation 
(11) are: 

C„,a(t) =C(0)[(1- {t/Tn,aT + ..] (16) 

independent of n for i << t^ q, while for n 7^ 1 and t » r„^Q, 

CuAt) = C(0)(n - 1) (VTn,a)-"/^"-'^ (17) 

The ratio of the the two asymptotic exponents a and a/{n — 1) yields the 
value of the apparent order of the reaction n. 

For special values of the two parameters n and a, some other typical solutions 
are recovered 



a. If n = 1, et = 1, we have 

dc{t) 



dt 

which is a first order kinetic 
h. li n — 1, a ^ 1, we have 

dCa{t) 



Kic{t) c{t) = c(0) exp{-Kit) (18) 



df" 



^K^c^it)^ c^{t)^c{0)eM-Kan (19) 



which is a "WeibuU kinetics". If < a < 1, it is a "stretched exponential 
kinetic" . 
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c. li n ^ 1, a — 1, equations (11) gives 

c,{t) = c„(0)[(l + (n - I)cn{0r-'KJ]-T^^ = expJ-c(0)"-iX„t) (20) 
which is solution of (2). 



d. If n = 2, Q! = 1, we have 
dc{t) 



K2c\t) ^ l/c{t) - l/c(0) = K^t (21) 



dt 

This is the second order kinetic, 
e. If n = 2, Q! 7^ 1, we have 

= K,,^c^{t? - c^{t) = c(0)[(l + cmK2,at)T' (22) 
This is a generahzed second order kinetic. 
Cases (b) and (e) have been discussed in [4]. 

It is important to note that as soon as n 7^ 1, the time dependence of the 
kinetics depends on the initial concentration. 

We will call the kinetic giving rise to the concentration evolutionary law (11), 
the (n, a) kinetic: 

CnAt) = c(0)[(l + {n-l) (i/T„,„)"]-^ (23) 

rn,a= [c{Or-'K^,a]-'^" (24) 

is the characteristic time of the complex kinetic. It depends on the initial 
concentration and the two exponents n and a. For n ^ 1, Cn^ait) tends to a 
WeibuU exponential with ti^q = [ii'i.„]^^/" . One can define a "half-reaction 
time" Ti/2 which is the time necessary to transform half of the relevant quantity 
by solving the equation 

(1 + (n - 1)(TV2/Tn,a)") -'^^""^'^ = 1/2 (25) 

which gives using the definition of ln„(a;) (4) : 

Ti/2 =T„,«(ln„2)V« (26) 

Kinetics are "memoryless" only when n = Q; = l. Ifn=l, kinetics are 
"memoryless" in the fractal time tf = since with the change of variable, 
the WeibuU function reduces to a memoryless exponential. 
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One can introduce in this problem a "response function" as it is done in Weron 
et al.[4] for n = 1 and n — 2. We have more generallyfor any real n: 

V / ' n,a ' n,a 

this function has the two asymptotic behaviors: 

for t^O f{t)0i{t/Tn,aT~^ (28) 
for t^OO (X (VTn.a)-^'^/"-^)-^ Tl^l 

For n — 1 and n = 2, they coincide with those of [4]. 
2 Arrhenius law and exponential conspiracy 

The results of the previous section can be understood physically as a con- 
sequence of what is has been called "exponential conspiracy", an expression 
coined by Boucheau [21] and proposed as exercise in textbooks on probabil- 
ity theory (for example [22]). It is generally accepted that the temperature 
dependence of the reaction rate K has an Arrhenius form which we will write: 

K ^ Koex.p(±E/kT) (29) 

Kq is the pre-exponential factor and E the relevant energy (in thermodynam- 
ics systems this energy is the Gibbs free energy which depends on the enthalpy 
(heat of reaction) and the entropy : G — H — TS) . The sign -|- corresponds to 
an "exothermic" reaction (i.e. the energy corresponds to an attraction energy 
and the rate decreases with the temperature). This is the case for instance in 
physisorption, when the overall adsorption enthalpy resulting from adsorption 
and desorption is positive. The sign - corresponds to an " endothermic" reac- 
tion and E is an activation energy barrier to be overcome. In that case the 
rate increases with temperature. We have written the two terms "exothermic" 
and "endothermic", because due to the variation of entropy with T, paradox- 
ically in some complex systems, an endothermic reaction can occur without 
activation energy (see for instance [23]). 

In disordered systems frozen out of equilibrium, the exponent factor E/kT 
varies due to fluctuations of local energies and local temperatures. The dis- 
tribution of energies depends on what has been called by Frauenfelder [3] the 
"energy landscape" , a concept taken from the theory of glasses, and variations 
of the inverse of local temperature (1/T) have been used to introduce, what 
has been called "super-statistics" by Beck and Cohen [24]. 

Here we will assume, as in the theory of heterogeneous catalysis, in the theory 
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of glass and in the theory of adsorption [25], that the probabihty distribution 
of the energy E varies for large values as exp(— i?/i?o)- This means that the 
large energy values are statistically exponentially small. The reference energy 
Eo is linked to the width of the energy density distribution fsiE) [25]. With 
this assumption, using the basic probability theory relation 

fE{E)dE = fK{K)dK (30) 

it is straightforward to show, that the distribution of the rate coefficient K, 
has the asymptotic form 

fK{K) = ii^K/Ko)-^^^' with ii = kT/Eo (31) 

Therefore, in the "exothermic" case (sign + in (18)), 

K ^oo fK{K)r^ K-^-" (32) 

while in the " endothermic" case (sign - in (18)) 

K ^ fK{K) ~ (33) 

In the first case, the distribution fxiK) is a Pareto distribution and is the 
simplest density distribution belonging to the domain of attraction of the sta- 
ble Levy distributions [26] [27]. If we assume that the variations of K induced 
by the fluctuations of E are represented by a Levy distribution L^{X) , one 
can, using the well-known relation (the Laplace transform of a one-sided Levy 
distribution is a stretched exponential) : 

roo 

/ exp{-XKt)C^, {X)dX^exp{-{Kt)^) (34) 

t/ 

obtain the generalized flrst order (WeibuU) kinetic (case 2 with = a) as a 
compounded exponential first order kinetic . In the second case, if we use the 
Gamma density distribution which has the power law asymptotic behavior 
(liX^"^) for small values of A, 

^'^^'^^ = 1^^^'^^'"'"''^^"^^^ ^^^^ 
we obtain equation (3) with = l/(n — 1) 

r eM-HKt))9,WdX = (1 + -{Kt))-'' (36) 
Jo /i 

If we use the WeibuU distribution, we can then write [28] 

r eM-HKtr)9,WdX = (1 + -{Ktry (37) 
Jo 11 
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or using a result established by Weron and collaborators in the stochastic 
theory of relaxation [6] : 

POD 

/ exp{-XKt) ML„^^{\)dX = (1 + {1/ ^){KtY)-^' (38) 

J 

where MLa,n{\) is a generalized Mittag-Lefller distribution: 

This last result is more difficult to interpret. It can be understood [6] [9], as 
the result of the random character of the number of active centers, geometric 
frustrations and dynamic constraints or as a consequence of the interplay of 
"exothermic" and " endothermic" processes in the kinetics of complex mate- 
rials. In conclusion, the (n, a) kinetic equation, can be obtained formally by 
introducing an adequate distribution for the exponent of the Arrehnius law as 
conjectured by Praunf elder [3]. 



3 Probabilistic interpretation of the (n, a) kinetic equation 



A comparison with the stochastic theory of relaxation [6] [4] [9] [10] is of interest 
to understand the physical meaning of equation (23). We first note that 

CnAt) = C(0)[(1 + (n - 1) (VTn,a)1"^ (23) 

is related to the Buyyxh distribution function {Ba,b,c{x) = 1 — (1 + ax'^)~'^ 
, a; > 0) [28] , named by reference to the number it occupies in the main Table 
of Burr's original paper [29]). If we introduce an effective random reaction 
waiting time 9 , the quantity c„,a(t)/c(0) can be viewed as the probability 
that the reactant has not yet reacted at time t: 

Cn,ait)/ciO) = Pr(6' >t)^l- f fnMde (40) 

J 

where 

fuAO) = + - (41) 

This distribution belongs to the domain of attraction of the Levy distribution 
with a tail exponent = a/{n — 1) and therefore generalizes the Pareto or 
Zipf-Mandelbrot distributions used in fractal reactions kinetics of previous 
works [5]. If /i < 1, an expectation value of 9 cannot be defined and an escort 
probability function [30] has to be used to determine r^^^ from the knowledge 

0ifn,a(9)- 
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If « = 1, the probability density function (41) reduces to the Tsalhs general- 
ized Pareto distribution [8] if (n — 1) = (g — l)/(2 — q) i.e. n — 1/(2 — g) 



m = —(1 + {q- l)/(2 - q)—)-- 

'n,a 'n,a 

which maximizes the Tsallis entropy of the random variable 6. 

The relation between the (n, a) kinetics and the nonextensivity of the en- 
tropy and the formal relation of the reaction order n with the Tsallis entropy 
index q is worth further investigations. The characteristics of the complex 
systems studied in the present work (see the introduction) are similar to the 
ones (frozen non-equilibrium states with memory effects and long range cor- 
relations) of what has been called nonextensive systems by the Tsalhs school 
[8]. 



4 Application to biotechnology and biophysics 



In this last section, we give some examples of problems in the field of biothech- 
nology and biophysics, where we think the application of the {n,a) kinetics 
can open new paths to understand anomalous kinetics from the point of view 
of the theory of complex systems. 



4-1 Sorption of dissolved contaminants in liquid phase 



The sorption (adsorption, chemisorption, biosorption) of pollutants from aque- 
ous solutions plays a significant role in water pollution control. It is therefore 
important to be able to predict the rate at which contamination is removed 
from aqueous solutions and how this rate depend on physical, chemical, bio- 
logical and environment variables in order to design an appropriate treatment 
plant. Sorption of dissolved contaminant is a complex phenomena caused by 
several mechanisms including London-van der Waals forces. Coulomb forces, 
hydrogen bonding, hgand exchange fluctuations, chemisorption, dipole-dipole 
forces and hydrophobic forces and biosorption for biological materials. There- 
fore these systems can be considered to belong to the class of complex systems 
[31]. The quantity adsorbed at time t, q^ is defined as 

where Cq is the initial concentration of the solution, Ct, the concentration at 
time t, V, the volume of the solution and W, the weight of the adsorbent. 
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In this context, the kinetic equations are determined with reference to the 
quantity of dissolved contaminant necessary to reach equilibrium(for ex. [23]). 
We have then 

qt = qe{l-exp{-K,t) (43) 

where Qe is the mass of solute adsorbed at equilibrium, qt is the mass of solute 
adsorbed at time t and Ki is the rate coefficient. Equation (43) is called pseudo 
first- order equation by contrast to the simple exponential first-order equation 
(18). 

In the same way, one can define a pseudo second-order reaction: 

^K2t + - (44) 

Qe - qt Qe 

In agreement with the ideas developed in sections 1 and 2, we can introduce 
the pseudo-(n,Q;) equation 

qt{a, n) = qe[l - (1 + qr\n - l)X„,«t")-^] (45) 

which reduces to (43) for n = 1, a = 1 and to (44) for n = 2, a = 1. We can 
write (45) more compactly using the definition (4) of the deformed exponential 
exp„(x), 

'^n,a(^) = ?t(",'T')M = l-exp„((t/r„,a)") with r„,a = (^e'^^n.a)"^^" 

(46) 

The definition of the deformed logarithm (4) associated with exp^(x) allows 
us to write the following relation 

R{t) = Log(^ ) = «Log(t) - aLog(r„,J (47) 

n — 1 

which can be used to make a linear fit of the data (rn,a(t) = qt{a,n)/qe) and 
obtain the values of a and Tn,a- The value of n to be chosen is the one which 
can give the better fit in the Log-Log plot. 

In two different collaborations we have analyzed the kinetics of various pol- 
lutants (phenol, tannic acid, gallic acid, melanoi'dine on activated carbon [32] 
and various dyes poUutants from the textile industry on biological materials 
(algaes and agaves) [33]. It appears that, quite generally the data can be fitted 
quite well to pseudo-(n,a) kinetic. For t << r, the concentration qe{n, a) does 
not depend on n (cf eq.l6) and can be fitted to a n- independent power law 
qe{n,a) cc t"' . The value of n (i.e. the overall order of the reaction) has to 
be determined from the large time (near saturation) behavior of the kinetics. 
Contrary to the usual practice in that field, we found that the exponent a, 
which is implicitly equal to 1 in the traditional analysis of kinetic data [23] 
in terms of a pseudo-first or -second order reaction, is more appropriate to 
characterize the kinetics of sorption of dissolved contaminants in liquid phase. 
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It yields a better fit and moreover it is related to the system energy land- 
scape which varies with external parameters. In Fig.l, we show an example 
(adsorption kinetic on activated carbon of melanoidine, a dye formed in the 
crystaUization process of sacharose [32]) where the best fit is obtained with 
a pseudo-(1.5,Q;) kinetic equation. One cannot fit properly the data with a 
quasi- first or second-order kinetic where a = 1. This is understandable since 
one observes that for small t , qe{n,a) oc t^-^^. The dependence of the two 
quantities a and n on the physical, chemical environmental and biological pa- 
rameters of the couple adsorbent-pollutant (pH, T, clustering, ligand field and 
architecture of large biomolecules...) is the subject of current studies [32][33]. 

4-2 Kinetics in photosynthesis processes 

In this subsection, we want to show how the fractal (n, a) rate equation (41) 
can be used in situation where a two-steps first order equation has been used 
to fit the kinetics of photosynthesis processes. The example chosen is the kinet- 
ics of the conversion of protochlorophyllide into chlorophyllide. The method 
used to follows the kinetics in that problem is the observation of spectral 
changes recorded by the technique of spectrofluorometry under short-time il- 
lumination. For instance in [34] the authors have observed the transformation 
of a 647 nm pigment (species a) by 630 nm photons. Contrary to the trans- 
formation induced by 647 nm photons, where the kinetics is first order, the 
transformation under 630 nm irradiation follows an unusual kinetic. It has 
been assumed that this particular kinetics is due to the presence of an other 
protochlorophyllide species (named b), i.e. a pigment with another associa- 
tion with the hpoproteins. We refer to the specialized literature for details 
[34] [35] [34] . The two-steps model, often used when the rate equation cannot 
be fitted to a first or second-order kinetics, lead Boardman, in this particular 
problem [35] , to the following rate equation (percentage of phototransformed 
quantity) , 

T{%) = 100 -Aex.p{-Kit) - (100 - A) ex.p{-K2t) (48) 

A is the proportion in % of the complex which is transformed. This formula 
can be fitted to the experimental data up to 85% [34]. What differs in 647 
pigment and 630 nm pigment is the link with the lipoproteins. The pigment- 
protein links are most probably fluctuating locally due to the complexity of 
the organization of the molecules inside the prolamellar body. If instead of 
two, a distribution of "species" is present, it is more appropriate to use, a 
rate equation deriving from a distribution of exponentials as the (n, a) rate 
equation. We have verified that in this particular case, (42) can flt perfectly 
the experimental curve, also for transformation larger than 85% (Fig.l). The 
two asymptotic behaviors appears to be power-law, a behavior a simple two- 
step mechanism cannot account for. The best fit obtained with a nonlinear 
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method is given by 

T(%) = 100(1 - (1 + (1.9 - 1) (t/0.07)°-96)-i/(^-^-i))) (49) 

This correspond to a fractal-time exponent a = 0.96 and a reaction order pa- 
rameter n = 1.9. The characteristic time is Tq, = 0.07 and from the inspection 
of the experimental curve (Fig.l), one can see that it is very close to Ti/2- The 
linear fit (r(t) = T/lOO ) gives the same result (Fig. 2): 

(1 - r(t)V-'^-^ - 1 

R{t) = Log{- — ) = 0.96Log(t) - 0.96Log(0.07) (50) 

1.9 — 1 

with a regression coefficient of 0.9997. 

The temperature dependence of the transformation rate under 633 nm photons 
[36] indicates that both endo- and exothermic effects are competing. As we 
have suggested at the end of section 2, in this particular situation, fluctuations 
of the rate coefficient in the exponent of the Arrhenius law (29) can give rise 
to {n, a) kinetics. 



4.3 Complexity of DNA 



One method widely used to study the complexity of DNA is the so-called Cot 
method. The method splits the double strands of DNA into single strands by 
raising the temperature or by other denaturing process. One then studies the 
kinetics of the reassociation of dissociated single strands. Since it involves two 
single strands, the renaturing into the original form is assumed to follow a 
(2,1) kinetics f = c (t)/c (0) = [(1 + c(0)(/^t)]^^. This equation is the basis 
of the Cot analysis of the rate of renaturation of sequence heterogeneity (or 
complexity) of DNA, The quantity c(0) is the initial concentration of DNA, 
/= c {t)/c (0) the fraction of single-stranded molecules which decreases with 
time and K the rate constant for the reassociation of complementary strands. 
The value of c(0)t when / = 0.5 is known as c(0)ti/2. 

The rate coefficient K is characteristic of a particular DNA and is related to 
its complexity in terms of sequence composition. The quantity c(0)ii/2 is the 
reciprocal of K and can therefore be used as a measure of sequence complexity. 
The higher the value of c(0)ti/2, the more complex is the DNA. 

The Cot method which was developed in the 1960 's and widely used in the 
70' was then nearly abandoned. It made a comeback recently [37] as a much 
cheaper method because of its ability to concentrate on the low copy sequences, 
the highly repeated sequences being irrelevant as far as the genetic information 
is concerned. 
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Due to the complexity of the DNA stnictTirc, it would be surprising if there 
would be no short time memory effects in the reassociation process, since 
it involves quite complex biomolecules. Indeed, fingerprints of fractality and 
nonextensivity in DNA fragment distribution has been reported recently [38]. 
The ideas and formalism developed in this paper might be of interest also in 
that field of primary importance. 



5 Discussion 



Using ideas and theoretical tools borrowed from recent works on the theory of 
relaxation, we have derived a universal function for the kinetics of complex sys- 
tems characterized by stretched exponential and/or power-law behaviors This 
kinetic function unifies and generalizes previous theoretical attempts to de- 
scribe what has been called "fractal kinetic". The concentration evolutionary 
equation (12) is formally similar to the Burrx/j relaxation function obtained 
in the theory of relaxation, with two exponents a and n. The first one is due 
to memory effects and short-range correlations and the second one finds its 
origin in the long-range correlations and geometrical frustrations which give 
rise to ageing behavior. As in the theory of relaxation, these effects can be 
formally handled by introducing adequate probability distributions for the 
rate coefficient. We have shown that the distribution of rate coefficients is the 
consequence of local variations of the free energy appearing in the exponent 
of the Arrhenius formula. The scaling (power-law) behavior of the kinetic is 
therefore an other example of what has been called "exponential conspiracy" 
[21]. The two macroscopic observables n and a are formally related to the 
energy landscape of the complex system which varies if physical, chemical or 
biological external factors are modified. 

The fractal (n, a) kinetic has been applied to a few problems of fundamental 
and practical importance [32] [33] [39] , examples of which have been presented 
in section 4. 

In references [9] we have shown how a universal relaxation function can be 
derived if we use distributions of macroscopic waiting times maximizing the 
nonextensive Tsallis entropy. Similar conclusions can be drawn in the present 
problem, if we introduce local reaction waiting time in a probabilistic deriva- 
tion of the universal kinetic function. 

The relation between the (n, a) kinetic and nonextensive thermostatistics will 
be the subject of further studies. 
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6 Captions of Figures 



Fig.l (n, a) kinetic of adsorption of melanoi'dine in aqueous solution on acti- 
vated carbon with n — 1.5 and a — 0.56. 

Fig. 2 (n, a) kinetic (eq.49) of the conversion of protochlorophylhde into chloro- 
phyUide (transformation of 647 nm pigment by 630 nm photons) . The smaller 
points for t>0.4 sec. are results of the two first order model (eq.48). 

Fig. 3 Log-Log plot (eq.47) applied to data of Fig. 2. 
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Fig.2 Kinetic ( eq.49) of the conversion of protochlorophyllide into chlorophyllide . 
The smaller points for t>0.4 sec. are results of the two first order model ( eq.48). 



